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ABSTRACT 

The  subject  of  this  paper  is  an  investigation  of  the 
significant  factors  that  affect  the  performance  of  a  sub- 
cavitating  pitot  type  intake  for  a  waterjet  propulsion 
system  for  hydrofoil  craft.  A  computer  model  for  design 
of  an  inlet  nacelle  is  developed  which  accepts  values  of 
two  parameters  and  several  system  constants  from  a  control 
program  and  designs  a  nacelle  to  meet  those  requirements 
subject  to  limitations  imposed  by  cavitation  and  geometry. 
The  parameters  are  inlet  velocity  ratio  and  inlet  diameter 
to  maximum  diameter  ratio.  A  body  of  revolution  config- 
uration with  profile  shape  based  on  certain  desirable 
pressure  distributions  is  assumed.  Extensive  use  is  made 
of  empirical  data  from  the  available  unclassified  literature. 

The  design  conditions  are  those  existing  at  maximum 
speed,  called  cruise  speed,  and  the  minimum  foil-borne 
speed,  called  take-off  speed. 

The  subroutine  will  accept  up  to  three  additional  speeds 
at  which  a  performance  prediction  will  be  made  after  design 
is  completed  from  th  first  two. 

Pressure  recovery  characteristics  and  cavitation  limits 
for  some  typical  designs  are  displayed  as  graphs. 

The  subroutine  listing  and  flow  diagram  is  included  in 
Appendix  C.   In  addition,  an  alternate  design  method  which 
accepts  only  the  inlet  velocity  ratio  as  a  parameter  is 
included  in  Appendix  D. 
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Chapter  1. 
PROBLEM  STATEMENT 

Ram  (pitot)  type  inlets  for  water jet  propulsion  systems 
in  hydrofoil  craft  are  responsible  for  a  significant  part 
of  both  the  external  drag  and  internal  flow  losses  and 
have  a  strong  influence  on  the  propulsive  efficiency  of 
the  system. 

The  objective  of  this  investigation  is  to  identify  and 
quantify  the  significant  characteristics  of  a  "well  de- 
signed" inlet  in  a  manner  that  will  allow  their  use  at  the 
preliminary  design  stage  to  assist  in  rapidly  determining 
the  feasibility  and  relative  merit  of  various  propulsion 
system  configurations  in  a  craft  of  a  given  size  with 
given  speed  and  range  requirements. 

The  term,  "well  designed",  as  used  with  reference  to 
the  inlet  and  nacelle  will  mean  that  inlet  configuration 
which  best  suits  the  propulsion  system  as  a  whole  when 
optimized  with  respect  to  some  chosen  criterion. 

A  computer  subroutine  is  to  be  developed  for  use  with 
an  optimization  program  concurrently  being  developed  by 
R.  P.  Gill  (ref.  1)  which  calls  this  subroutine  and  others 
written  by  Gill  and  R.  C.  Percival  (ref.  2).  The  optimizat- 
ion criterion  chosen  for  this  particular  program  is  least 
propulsion  system  weight  including  fuel,  but  the  inlet 
subroutine,  called  SUBROUTINE  NACEL,  does  not  depend  on 
that  particular  criterion  and  could  be  used  for  others. 

Engineering  approximations  and  empirical  fits  are 
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used  wherever  they  are  suitable  and  consistant  with  the 
purpose  of  the  study. 

There  are  at  least  three  broad  categories  of  nacelles 
which  could  be  used  for  ram  inlets:  subcavitating,  super- 
cavitating,  and  base-vented.  Only  subcavitating  nacelles 
are  considered  in  this  study.  Persons  who  are  interested 
in  supercavitating  and  base-vented  nacelles  may  wish  to 
read  references  (3)  and  (4)  for  further  discussion  on 
those  types. 

Subsonic  aerodynamics  provides  much  of  the  data 
utilized  in  this  study,  but  hydrodynamic  design  has  the 
additional  complication  of  cavitation  avoidance.  For 
the  purpose  of  this  paper  cavitation  is  assumed  to  start 
when  the  local  static  pressure  decreases  to  the  vapor 
pressure  of  sea  water.  Incipient  cavitation  is  imposed 
as  an  absolute  limit  on  the  decrease  in  static  pressure 
at  the  design  conditions  even  though  a  certain  amount  of 
cavitation  may  be  allowable  in  a  real  life  situation. 


-12- 


Chapter  2. 
DESIGN  CONDITIONS 

There  are  three  possible  flow  conditions  at  the  inlet 
and  Bernoulli's  equation  indicates  the  changes  in  local 
static  pressure  coefficient  associated  with  each  of  them. 

When  Vi/Vo-^1,  the  stream  lines  are  diverging  (Fig.l) 
and  since  there  are  no  losses  (to  the  inlet  plane,  at 
least),  local  velocities  inside  are  lower  than  free  stream 
velocity  causing  the  pressure  coefficient  to  be  positive 
and  no  cavitation  danger  exists.  Local  velocities  on  the 
exterior  of  the  forebody  are  higher  than  in  the  free  stream 
and  the  pressure  coefficients  are  negative  (Fig. 2),  so 
cavitation  danger  exists  and  increases  as  Vi/Vo  decreases. 
Thus,  for  any  particular  nacelle  and  free  stream  velocity 
combination,  there  is  a  minimum  velocity  ratio  below  which 
cavitation  occurs  in  the  external  flow. 

When  Vi/Vo^-1,  the  stream  lines  are  converging 

(Fig.  1).  Local  velocities  in  the  external  flow  will 

still  be  somewhat  greater  than  in  the  free  stream  but 

lower  than  in  the  diverging  flow  case,  so  if  the  external 

flow  was  cavitation  free  when  Vi/Vo  **1,  it  will  also  be 

cavitation  free  when  Vi/Vo  ^»  1  at  the  same  Vo.  Local 

velocities  in  the  internal  flow  are  now  greater  than  the 

free  stream  with  a  corresponding  decrease  in  static  pressure 

coefficient,  and  cavitation  inside  the  lip  is  a  definite 

possibility.  There  is  therefore  a  maximum  velocity  ratio 

for  each  nacelle  and  Vo  above  which  cavitation  will  occur 
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Figure  1.  Entering  stream  tube  illustrating  two  flow 
conditions  into  the  nacelle. 
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Figure  2.  Minimum  C  in  the  external  flow  vs.  V./VQ  for 


several  of  the  Lockheed  nacelles 
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in  the  internal  flow. 

When  Vi/Vo  =  1,  the  flow  is  in  an  intermediate  state 
between  the  other  two  conditions  and  no  new  limitations 
are  imposed.  This  condition  is  of  no  interest,  then,  from 
the  design  viewpoint. 

Changes  in  local  velocity  due  to  the  geometric  con- 
figuration of  the  nacelle,  lip  shape,  bluntness  of  the 
forebody,  etc.,  are  perturbations  on  the  flow  conditions 
described  above  and  are  assumed  to  be  minimized  by  the 
detailed  design  of  the  nacelle.  No  attempt  is  made  here 
to  define  that  detailed  design  except  to  assume  that  it  is 
at  least  as  good  as  the  best  of  the  empirical  examples 
and  data  available  to  this  investigation. 

A  very  frequently  encountered  set  of  requirements, 
as  indicated  by  the  literature,  includes  the  following 
situation: 

Cruise  speed  is  approximately  twice  take-off  speed. 
Cruise  flow  rate  is  approximately  equal  to  take-off 
flow  rate. 
Let  the  subscripts  1  and  2  designate  cruise  and  take-off 
conditions,  respectively.  If  the  inlet  area  is  a  constant, 

Q^Vij  Ai  =  Q2=Vi2  Ai 
Vi^Vig 
but  702=702-2,  so  Vi/Voi=^Vi/Vo2 . 
If  Vi/Vo^O.?  to  0.£  as  is  frequently  the  case, 
(refs.  4,5,6,)  then  Vi/Vo2=1.4  to  1.6,  but  except  when 

the  inlet  is  very  deeply  submerged,  these  high  velocity 
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ratios  will  cause  cavitation  inside  the  inlet,  incurring 
high  losses  and  perhaps  choking  the  inlet. 

One  way  of  decreasing  this  effect  is  to  decrease  the 
inlet  velocity  ratio  at  cruise  by  increasing  the  area. 
This  approach  very  quickly  reaches  the  minimum  velocity 
ratio  below  which  cavitation  occurs  in  the  external  flow 
with  its  accompanying  high  drag  and  possible  material 
damage . 

As  speed  increases,  the  interval  between  these  two 
limits  becomes  smaller,  so  the  problem  grows  increasingly 
difficult  as  attempts  are  made  to  achieve  higher  speeds. 

A  solution  to  the  problem  may  be  found  in  the  use  of 
auxiliary  inlets  to  increase  the  effective  inlet  area  when 
high  velocity  ratios  are  required. 

If  the  existance  of  an  auxiliary  inlet  is  postulated, 
the  velocity  ratio  at  take-off  is  no  longer  coupled  to 
the  velocity  ratio  at  cruise  and  both  become  design  con- 
ditions as  the  total  inlet  area  for  each  must  be  determined 
independently. 

An  additional  limitation  is  imposed  by  the  elbow 

where  the  internal  flow  leaves  the  nacelle  and  enters  the 

strut.  This  elbow  will  normally  contain  turning  vanes  to 

help  minimize  losses.  Since  the  local  velocities  on  the 

turning  vane  surface  will  be  greater  than  the  average 

velocity  in  the  flow  entering  the  elbow,  a  danger  of 

cavitation  exists.  This  turning  vane  cavitation  is  avoided 

by  diffusing  the  flow  in  the  nacelle  to  reduce  the 
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velocities  and  increase  the  static  pressure.  Therefore, 
a  minimum  diffusion  ratio  dictated  by  a  maximum  nacelle 
exit  velocity,  above  which  cavitation  occurs  in  the  elbow, 
is  a  third  design  limit.  This  limit  must  be  checked  at 
both  cruise  and  take-off  because  it  is  a  function  of  both 
the  flow  rate  and  the  losses. 

In  summary,  then,  both  cruise  and  take-off  are  design 
conditions  and  cavitation  avoidance  establishes  absolute 
limits  on  the  minimum  velocity  ratio  at  cruise,  the  maximum 
velocity  ratios  at  both  cruise  and  take-off,  and  the 
minimum  diffusion  ratio. 
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Chapter  3. 
NACELLE  SHAPE 

The  principles  of  design  used  in  jet  engine  nacelles 
and  air  intakes,  already  established  in  the  field  of 
aerodynamics,  can  be  used  for  guidance  in  the  selection 
of  waterjet  inlet  nacelles  (ref .  3  and  4) • 

Kiichemann  and  Weber  (ref.  7)  have  developed  the  de- 
finition of  a  theoretical  optimum  nacelle  contour  by  con- 
sidering the  momentum  balance  between  the  free  stream  and 
the  entry  plane  of  a  semi-infinite  intake.  The  resultant 
force  on  the  nacelle  is  given  by 

F  ■  i  V2  A.  (1  -  V./V  )2, 
N      o  i       1'  o   ' 

but  since  there  can  be  no  resultant  force  on  a  semi-infinite 
body  in  potential  flow,  this  force  must  be  equal  and 
opposite  to  the  integral  of  the  pressure  over  the  forebody 
area  of  the  nacelle. 

V  -^o/   CP  dA  =  -i/^oy  U  -  <W21 dA 

If  the  pressure  coefficient  is  a  constant  over  the 
entire  forebody  then  the  local  velocity,  V,,  is  a  constant 
and  the  projected  frontal  area,  Am-A.,  required  to  balance 
F„  is  defined  by  the  above  equation.  When  the  local 
velocity  is  limited  to  some  maximum  value  to  satisfy  Mach 
number  requirements  or,  in  the  case  of  waterjet  nacelles, 
cavitation  requirements,  then  the  minimum  projected  frontal 
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area  and  V  will  satisfy  the  relationship 

max*  '  r 

*£i1+  (1  -  Vi/Vp)2 
Ai       (Vnax/VoF  -  1 

The  optimum  intake  contour  of  Kiichemann  and  Weber  is 
one  which  at  its  design  V./V  achieves  the  goal  of  keeping 
Vmax  constant  and  at  its  specified  value  over  the  entire 
forebody  and  the  equals  sign  applies. 

A  /A.  can  be  represented  by  the  square  root  of  its 
inverse,  D^/D  ,  in  the  case  of  circular  intakes. 

(1  -  Vi/V0)2 


Pi/Dm  - 


1  * 


v  max'  o'   -1- 

It  can  be  seen  that  a  maximum  value  of  D./D  is  defined 

l'  m 

by  a  minimum  value  of  V^/V  for  any  given  VQ.   One  of  the 
design  limits  chosen  in  the  previous  chapter  was  avoidance 
of  cavitation  in  the  external  flow  at  the  relatively  low 
velocity  ratio  at  cruise  speed.  This  can  now  be  translated 
into  a  maximum  value  of  D^/Dm. 

It  is  not  practical  to  build  such  perfect  intakes 
although  intakes  with  reasonably  uniform  pressure  dis- 
tributions at  certain  velocity  ratios  have  been  developed 
(ref.  3,  U,   7).  In  general,  the  velocity  distribution 
will  have  peaks. 

The  generation  of  practical  contours  with  acceptiable 
pressure  distributions  has  been  largly  a  matter  of 
experimental  trial  and  error.  Kiichemann  and  Weber  de- 
scribe three  sets  of  geometrically  similar  shapes  which 
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they  designate  as  Classes  A,  B,  and  C.  The  geometric 
description  of  these  are  given  in  Table  1  and  Figure  3 • 

References  3  and  4  describe  two  sets  of  intakes 
whose  pressure  distributions  were  generated  by  numerical 
solutions  to  the  potential  flow  problem.  The  Lockheed  in- 
takes were  chosen  for  input  data  for  this  study  chiefly 
because  the  data  covered  a  wide  range  of  geometric  variation 
and  data  on  the  interal  flow  distribution  was  available 
for  each  of  the  individual  inlets.  The  major  geometric 
parameters  and  peak  pressure  coefficients  for  each  of  the 
nine  inlets  is  given  in  Tables  2  through  11.   The  external 
contour  of  these  nacelles  is  very  close  to,  but  not  quite, 
elliptical.  The  internal  lip  contour  is  elliptical.  The 
Kiichemann  and  Weber  intakes  are  elliptical  both  internally 
and  externally.  Brown  and  Traksel  state  that  this 
variation  has  an  influence  that  is  secondary  to  the  length 
to  diameter  ratio,  L/D  ,  and  D^/Dm.  Figure  4  shows  the 
similarity  of  the  major  parameters.  The  dashed  line  was 
tabulated  as  the  basic  geometric  description  of  nacelle 
forebodies  for  SUBROUTINE  NACEL.   The  three  data  points 
to  the  left  represent  intakes  which  had  lower  pressure 
coefficients  than  the  others,  so  the  dashed  line  represents 
the  best  of  the  available  data  points. 

To  obtain  a  functional  relationship  for  the  variation 
of  pressure  coefficient  with  geometery,  the  peak  pressure 
coefficient  at  each  velocity  ratio  and  at  each  angle  of 
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attack  was  plotted  against  the  value  of  Dj/D  and  also 

against  L/Dm*  ^/®m   was  cnosen  as  tne  appropriate  parameter 

to  use  because  it  gave  a  smoother  distribution  of  data 

points.   The  desired  relationship  was  then  obtained  by 

drawing  a  smooth  envelope  that  enclosed  all  points  and 

passed  through  the  points  of  highest  C  .   This  envelope 

is  considered  to  represent  the  "well  designed"  inlet 

mentioned  in  the  problem  statement  based  on  the  assumption 

that  if  a  nacelle  has  already  been  designed  which  falls 

on  that  line  (the  discrete  data  points  through  which  the 

line  was  drawn),  then  a  nacelle  could  be  designed  for  any 

other  point  on  that  line.  The  plots  are  shown  in  Figures  5 

through  10.  Figures  11  through  16  show  similar  plots  for 

the  internal  flow  conditions. 

There  can  be  no  losses  in  the  stream  tube  from  out 

in  the  free  stream  up  to  the  inlet  plane,  but  from  the  inlet 

plane  on  into  the  ducting  of  the  nacelle  friction  and 

turbulence  take  a  toll  of  the  head  available  in  the  internal 

flow.  The  amount  of  head  loss  is  primarily  a  function  of 

velocity  ratio,  Reynolds  number  and  lip  shape.  Blackaby 

and  Watson  investigated  experimentally  several  lip  shapes 

including  circular  profiles,  elliptical  profiles  and  sharp 

profiles.   In  general,  the  elliptical  profile  caused  less 

losses  at  a  given  Vj/V  than  the  others.  The  particular 

results  for  the  lip  which  they  designated  1$E  is  given  as 

a  function  of  Reynolds  number  for  several  velocity  ratios 

in  reference  3,  Figure  74  of  Report  No.  5.  The  loss 
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coefficient,  defined  as  A  P/q  ,  approaches  a  constant  value 
for  each  velocity  ratio  for  Reynolds  numbers  greater  than 
&xl(A.  These  constant  values  as  read  from  the  citation 
above  are  given  in  Table  12.  The  Ram  Pressure  Recovery, 
defined  as 

RPR=  (P  -  p0)/(P0  -  p0)  =  1  -  AP/q0 
is  plotted  in  Figure  17,  and  points  from  that  plot  are 
tabulated  in  the  subroutine. 
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Chapter  4 
THE  DIFFUSER 

The  inlet  nacelle  modeled  in  this  study  is  assumed  to 
consist  of  a  streamlined  body  of  revolution  containing 
a  nose  intake  opening,  a  diffuser  and  an  auxiliary  inlet. 
The  forebody  and  intake  lip  serve  to  guide  the  flow  into 
the  internal  ducting.  That  internal  ducting  includes  the 
diffuser  which  has  at  least  two  primary  functions.   One 
is  to  slow  the  flow  velocity  to  below  the  critical  value 
for  cavitation  in  the  strut  elbow.  The  other  is  to  de- 
crease the  flow  velocity  further  to  reduce  ducting  losses 
if  it  benefits  the  system  as  a  whole. 

Three  main  profile  shapes  of  diffusers  with  circular 
cross  section  can  be  considered.  They  are  bell,  trumpet 
and  straight  walled  conical  diffusers.  To  avoid  confusion 
of  terms,  a  bell  diffuser  is  one  having  its  mean  diameter 
greater  than  the  mean  diameter  of  a  straight  walled  conical 
diffuser  with  the  same  inlet  and  exit  areas  while  a  trumpet 
diffuser  has  a  mean  diameter  less  than  that  of  the  cone. 
At  similar  inlet  flow  conditions,  the  bell  shape  results 
in  a  lower  total  head  loss  than  the  other  two  and  the 
trumpet  results  in  a  higher  loss  (ref .  4,  9) • 

In  general,  the  diffuser  with  the  least  loss  in  total 
head  is  preferred  for  this  application  if  other  con- 
siderations are  equal  because  the  difference  in  losses  must 
eventually  be  made  up  by  expenditure  of  prime  mover  power. 
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Current  knowledge  of  boundary  layer  growth  and  be- 
havior in  adverse  pressure  gradients  does  not  allow  accur- 
ate prediction  of  diffuser  performance  by  analytic  methods 
(ref.  3),  and  very  little  empirical  data  on  bell  diffuser 
performance  that  is  directly  applicable  was  found.  Ref- 
erence 9  contains  an  experimental  comparison  of  the  three 
shapes  which  indicates  that  conical  diffusers  have  pressure 
recovery  characteristics  very  close  to  the  bell  and  sig- 
nificantly better  than  the  trumpet.  Experimental  results 
contained  in  references  9  through  13  can  be  applied  to 
predict  the  losses  of  straight  walled  conical  diffusers 
and  can  be  considered  a  lower  bound  of  that  which  can  be 
expected  from  a  bell. 

The  total  pressure  loss  for  a  conical  diffuser  can 
be  predicted  by 

AP  =  (C  Kt  4  f  L/dm)  qlf   (ref.  13), 
where  the  second  term  is  a  ducting  loss  of  a  duct  the  same 
length  as  the  diffuser  and  of  a  diameter  equal  to  the  mean 
diameter  of  the  diffuser.   For  a  cone, 

dm  =  (D-l  4  D2)/2   . 
The  first  term  accounts  for  diffusion  losses  as  a  percent- 
age of  the  losses  in  a  sudden  expansion.   K^  is  the  sudden 
expansion  loss  factor  and  is  given  by 

Kt  =  (1  -  A-j/A,,)2   . 

The  proportionality  factor  C  is  a  function  of  ex- 
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pansion  angle,  0,  or  the  double  angle,  20,  and  is  determined 
experimentally  as  in  reference  13 • 

The  minimum  loss  for  any  given  diffusion  ratio  occurs 
when  the  value  of  20  is  in  the  vicinity  of  5°  to  10°.  Flow 
separation  starts  interfering  with  the  flow  and  partially 
stalling  the  exit  when  20  is  in  the  vicinity  of  20° 
(ref.  ll),  therefore  the  region  of  interest  is  for  20  less 
than  20°.   In  this  region  the  curve  of  C  vs.  20  can  be 
approximated  by 

C  =  3.190  x  10"3  (0)2  +  g.452  x  10 "^  (0) 
where  0  is  in  degrees. 

There  appears  to  be  no  reason  in  this  application  to 
decrease  the  angle  below  the  value  of  minimum  loss  because 
the  diffuser  length  increases  and  friction  losses  start  to 
rise  accordingly.  For  the  computer  model  limits  of  3°  and 
10°  were  set  on  the  minimum  and  maximum  values  of  ©  which 
has  the  effect  of  establishing  a  maximum  and  minimum  length. 

The  desired  total  length  of  the  nacelle  is  generally 
established  by  external  hydrodynamic  flow,  structural  and 
other  considerations  that  do  not  include  the  diffuser  and 
other  internal  length  requirements.  If  the  internal  com- 
ponent lengths  are  less  than  the  desired  length,  they  will 
all  fit  within  the  envelope.  If  not,  additional  external 
drag  due  to  the  added  length  and  wetted  surface  is  generated 
and  the  power  required  to  overcome  that  drag  increment  can 

be  charged  to  the  diffuser  length  since  other  internal 
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lengths  are  fixed.  This  becomes  a  criterion  for  determin- 
ing the  .proper  diffuser  length.  The  sum  of  the  power  re- 
quired to  overcome  diffuser  losses  and  the  additional 
power  due  to  the  increment  in  drag  is  minimized  with  re- 
spect to  diffuser  length.  The  calculations  used  in  this 
determination  are  given  in  Appendix  B. 

A  limit  on  the  minimum  diffusion  ratio  is  given  by  the 
requirement  that  the  maximum  diffuser  exit  velocity  must 
be  less  than  the  critical  velocity  which  will  cause 
cavitation  in  the  elbow.  Experimental  results  in  reference 
11  give  exit  velocity  profiles  for  six  diffuser  examples, 
in  two  of  which  26  =  10°  and  in  four  of  which  26  =  12°. 
The  diffuser  area  ratio  was  2.0  in  all  cases,  entry  Mach 
number  was  0.26,  and  boundary  layer  thickness  ranged  from 
0.003  to  0.03  times  the  entry  radius.   Integration  of 
these  distributions  graphically  gives  a  range  of 
^max./^ave.  =  ^••^■5  to  1.50,  with  a  mean  value  of  1.33 
over  the  six  cases.  However,  since  boundary  layer  thick- 
ness has  poor  predictability,  the  average  exit  velocity  is 
not  known  but  it  will  be  somewhat  greater  than  the  flow 

rate  divided  by  the  exit  area.   If  V__„  is  calculated  by 

max . 

k  =  1.33  will  probably  give  a  low  value  and  k  =  1.50  may 
give  a  high  value  but  will  be  on  the  safe  side. 

The  upper  limit  on  Vfflax#  will  be  the  critical  value 
of  V2  determined  from  Appendix  A. 
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Chapter  5. 

THE  AUXILIARY  INLET 

It  was  indicated  in  chapter  2  that  an  auxiliary  inlet 
was  probably  required  to  obtain  the  flow  rates  necessary 
at  take-off  speed  and  remain  subcavitating  at  cruise. 

The  possible  configurations  for  providing  auxiliary 
inlets  are  many  and  varied.  Slots  in  the  nacelle  wall  with 
doors  that  open  inward,  leading  edge  openings  on  the  strut 
and  a  movable  nacelle  lip  are  all  possibilities,  to  mention 
only  a  few. 

The  model  for  this  study  was  assumed  to  be  of  the  first 
type.   Figure  3  shows  the  auxiliary  inlet  as  dashed  lines 
in  the  profile.   It  is  assumed  that  the  inlet  is  in  the 
form  of  a  series  of  openings  in  an  annular  ring  around 
the  nose  opening  with  some  fraction  of  that  ring  being 
structure  to  support  the  nose.  The  dimension  X  in  the 
figure  is^  the  width  of  that  annular  ring  in  the  direction 
normal  to  the  inlet  flow  through  the  opening. 

The  auxiliary  inlet  flow  is  assumed  to  enter  and 
combine  with  the  flow  from  the  nose  opening  prior  to  enter- 
ing the  diffuser  because  that  is  the  area  of  lowest  static 
pressure  and,  since  the  inlet  will  only  be  used  at  take-off 
when  the  flow  rate  is  high,  it  is  important  that  the  entire 
combined  flow  pass  through  the  diffuser  to  prevent 
cavitation  in  the  elbow. 

When  the  inlet  is  closed  the  length  of  ducting,  L  , 
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associated  with  it  is  considered  to  be  a  short  length  of 
pipe  and  the  loss  for  that  length  of  pipe  is  added  to  the 
other  system  losses. 

This  type  of  inlet  appears  to  be  in  much  the  same  flow 
conditions  as  the  "gliding  scoop"  described  in  reference  14. 
Since  the  auxiliary  inlet  is  open  only  when  the  system  de- 
mands a  high  flow  rate  its  inlet  velocity  ratio  is  estimat- 
ed to  be  about  1.0  or  slightly  higher.   In  this  condition 
the  lip  loss  coefficient  is  approximately  0.2  (ref.  14). 
This  value  is  the  one  used  in  the  subroutine  for  auxiliary 
inlet  losses. 
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Chapter  6. 
DESIGN  RATIONALE 

The  design  of  an  inlet  for  the  specific  conditions 
imposed  by  the  optimization  program  is  based  on  comparison 
of  the  characteristics  implied  by  the  system  parameter 
values  with  those  stored  in  the  data  tables  and  sizing  the 
individual  components  to  remain  cavitation  free  at  the 
design  conditions.   The  subroutine  listing  and  a  computer 
generated  flow  diagram  are  in  Appendix  C.  This  and  the 
following  chapter  will  attempt  to  clarify  the  reasoning 
used  in  its  derivation. 

When  the  program  calls  SUBROUTINE  NACEL,  all  physical 
constants,  craft  configuration,  depth  of  submergence, 
free  stream  velocity,  flow  rate  and  angle  of  attack  for 
each  design  point  are  known.   Inlet  velocity  at  cruise 
and  D./D  are  passed  as  system  parameters.  Actually, 
V./V  at  cruise  is  the  system  parameter,  but  it  is  cal- 
culated from  V.  and  V  rather  than  passed. 

Initially,  free  stream  static  pressure,  free  stream 
dynamic  pressure,  incipient  cavitation  number,  and  total 
pressure  for  both  cruise  and  take-off  speeds  are  calculated. 

CToi  =  (P0  "  Py'i/^o 

Po  =  P0  +  % 

The  length  to  diameter  ratio,  L/D  ,  of  the  nacelle  for 

'   '  m' 

-30- 


which  the  external  minimum  pressure  coefficient  is  equal 

to  the  negative  of  the  free  stream  incipient  cavitation 

number  at  cruise  is  determined  by  interpolation  in  the 

data  table  where 

L/D  =  f(C  ,  V./V  ,  *). 
'  m    v  p'   i'  o' 

The  inlet  to  maximum  diameter  ratio,  D./D  ,  correspond- 
ing to  that  L/Dm  is  found  in  the  tabulated  geometric 
characteristics.   The  value  obtained  defines  a  nacelle  for 
which  the  external  flow  on  the  surface  is  at  the  point  of 
incipient  cavitation  and  therefore  is  a  maximum  value. 

(0i/Vraax  =  f<L/Dra> 

If  the  Di/Dm  passed  by  the  calling  program  is  greater 

than  (D^/D  )    ,  that  particular  set  of  parameters  is 

rejected  by  setting  the  nacelle  weight  to  a  very  high 

number  and  returning  to  the  calling  program.   If  D^/D  is 

less  than  (^i/Dm)max»  then  it  is  cavitation  free  and  the 

design  proceeds. 

Inlet  area  is  determined  from  the  flow  rate  and  inlet 

velocity  at  cruise  which  then  allows  D.  and  D  to  be 

1      m 

determined. 

The  maximum  allowed  velocity  ratio  for  both  speeds  is 
determined  by  interpolation  in  the  data  table  again. 

(V./V  )    =  f (L/D  ,  <r  .,  <*) 

i'  o  max      '  m*   oi* 

This  limit  is  imposed  by  cavitation  on  the  inside  of  the  lip. 

If  at  cruise,  V./V  is  greater  than  (V./V  )    ,  the 
'   i'  o    to  l'  o  max* 

parameters  are  rejected  in  the  same  manner  as  before.   If 
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not,  the  take-off  condition  is  checked  to  determine  if 

an  auxiliary  inlet  is  required.  The  maximum  flow  rate 

allowed  by  (V./V  )    enters  through  the  nose  inlet,  the 
J        r  o'max  &  ' 

remainder  must  enter  through  the  auxiliary  inlet. 

The  pressure  recovery  of  the  lip  is  a  tabulated 

function  of  V./V  ,  so  the  total  pressure  inside  the  inlet  is 
i'  o*  r 

Pi  =  RPRlip  ^0  +  Po 
and  the  static  pressure  inside  is  given  by 

p.  =  P-  -  kfV*    . 
Pressure  recovery  of  the  auxiliary  inlet  is  known 
from  chapter  5  so  the  total  pressure  inside  the  auxiliary 
inlet  is 

Paux  =  RPRaux  %   *  P0  • 
The  static  pressure  at  a  point  is  unique  so  the  static 

pressure  inside  the  auxiliary  inlet  is  the  same  as  that 

previously  calculated  from  the  flow  inside  the  nose  inlet. 

Therefore,  the  velocity  through  the  auxiliary  inlet  is 

V    =  [2(Paux  -  Pi)//7]*   • 
aux    L   &ux  r  l  '  J 

The  area  of  the  auxiliary  inlet  is  found  by 

Aaux  =  Qaux/  Vaux  • 
The  total  pressure  of  the  combined  flow  is  estimated 

as  the  mass  weighted  average  of  the  combining  flows. 

Pc  =  <W  Paux  *-Qi  Pi)/  Qc 
where  Qc  =  Qaux  +  Q±      . 

The  average  velocity  of  the  combined  flow  is 
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calculated  from  the  dynamic  pressure 

vi  =  [  2(Pc  -  Pi'//>]  4  , 

and  the  pressure  recovery  coefficient  to  this  point  is 

RPRC  =  (Pc  -  pQ)/  qQ  . 

The  diffuser  entry  and  exit  diameters  are  set  by 

Dl  =  Di  .  D2  =0.9  Dm  . 

The  0.9  factor  is  somewhat  arbitrary.   Some  fraction 
of  the  nacelle  diameter  must  be  allowed  for  structural 
requirements  and  it  was  felt  that  ten  per  cent  would  be 
sufficient  in  most  cases. 

The  internal  component  length  requirements  are 
calculated  in  the  manner  of  Appendix  B  and  the  sum  of 
those  lengths  plus  the  lip  length  is  used  as  the  nacelle 
length  subject  to  a  minimum  value  of  5*5  Dm»   This  factor 
is  also  arbitrary  and  is  the  average  of  the  recommended 
values  of  LN/Dm  found  in  the  literature.  Hoerner  (ref.  15) 
shows  from  empirical  data  that  the  minimum  drag  nacelle 
on  a  wing  has  a  L^/B     value  of  approximately  3*0,  but  since 
the  internal  length  requirements  are  always  greater  than 
3.0  (see  Appendix  B),  it  was  felt  that  value  might  place 
an  unncesssarily  high  drag  penalty  on  the  diffuser. 

The  larger  of  the  two  flow  rates,  cruise  or  take-off, 
is  used  for  sizing  the  diffuser. 

The  diffuser  loss,  lip  loss,  and  pipe  loss,  if  any, 

are  calculated  for  both  cruise  and  take-off.  Their  sum  in 

each  case  is  used  in  the  equation  derived  in  Appendix  A 
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to  determine  the  critical  velocity  value  at  the  diffuser 
exit  for  cavitation  on  the  elbow  turning  vane. 

Vcrit  "  koi  *  X  "  Ploss/i/'Vo  J  *  V(1  +  Ivi)* 

The  maximum  diffuser  exit  velocity  as  determined  in 
chapter  4  is 

V2„,ax  =  X'5  Q/A2  ' 
These  two  velocities  are  compared  for  both  conditions. 

If  V     is  greater  than  V     for  either  condition,  the 
4~max  criT3 

subroutine  returns  to  the  calling  program  as  before  with 
ai  arbitralily  high  weight  value.   If  V     is  satisfactory, 
the  design  is  complete. 

External  drag  is  estimated  by 

D  =  CD  i/^Vo  A 
where  A  is  the  external  surface  area  of  the  nacelle.  A  is 

determined  by  dividing  the  nacelle  in  to  three  sections, 

the  forebody  from  the  nose  to  the  point  of  tangency  with 

a  line  parallel  to  the  centerline  at  maximum  diameter,  a 

cylinderical  midbody,  and  a  tail  section  which  is  assumed 

to  be  symmetrical  with  the  forebody.   The  forebody  and 

tail  can  be  approximated  as  paraboloids  of  revolution 

and  their  combined  area  is 

Ap  =  (W3)(Dm/2)2  [M  -  1/(1  -  M)] 

where  M  =    [  1  -   4(2L/Dm)2]  4  (ref.   4). 

For  the  cylinderical  midbody, 

Ac  =  T7(Dm/2)2  (LN  -  2L) 
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and  then  A  =  A  +  A„ 

P    c 

The  drag  coefficient  is  from  Hoerner  (ref.  15)  and 
is  given  by 

CD  =  0,  [  1  ♦  1.5(Dm/V3/2  +  7(Dm/LN>3]   • 
The  Schoenherr  friction  coefficient,  C~,  is  calculated 

by  a  function  subroutine,  CFS(R),  in  which  the  argument 

is  Reynolds  number. 

The  weight  estimate  depends  on  materials  and  details 
of  the  structure  which  were  not  the  concern  of  this  pro- 
ject.  Therefore,  only  a  rough  approximation  is  possible 
and  the  assumptions  are  certainly  open  to  question. 

Internal  ducting  weight  was  specified  to  be  propor- 
tional to  the  cross  sectional  area  and  the  length  such 
that     W  =  15.07  (Area)  (Length)   in  pounds 
where  area  is  in  square  feet  and  length  in  feet  (ref.  1). 

The  weight  of  the  outer  envelope  is  estimated  by 
assuming  that  it  is  proportional  to  the  surface  area, 
the  maximum  wall  thickness,  and  the  specific  weight. 
A  maximum  wall  thickness  of  0.05Dm  was  used  to  be  con- 
sistant  with  the  diffuser  exit  diameter  assumption. 
The  outer  envelope  of  the  nacelle,  of  course,  is  not 
solid  but  is  made  of  inner  and  outer  skins  plus  stiff- 
ening structure.   The  specific  weight  depends  on  the 
material  which  has  not  been  specified.  Two  extremes  of 
the  weight  spectrum  which  might  be  used  for  such  a 

structure  are  aluminum  and  steel.  Some  material  with 
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a  weight  between  these  two  is  more  likely  and  a  material 
density  of  one  half  that  of  steel  was  used  in  the  sub- 
routine.  The  actual  weight  in  water  may  very  well  be 
close  to  zero  if  the  voids  are  filled  with  oil  or  pitch. 

The  proportionality  constant,  a  measure  of  the  vol- 
ume of  the  material  in  the  wall  compared  to  the  volume 
of  the  wall  itself,  a  packing  factor,  will  be  a  small 
decimal  on  the  order  of  0.01  to  0.1.   0.03  was  used  in 
the  subroutine. 

The  resulting  weight  was  increased  by  ten  percent 
to  allow  for  the  auxiliary  inlet  mechanism. 

At  this  point  the  subroutine's  job  is  completed  and 
it  returns  to  the  calling  program  with  the  following 
quantities: 

a.  the  head  loss  summed  from  the  free  stream  to  the 
diffuser  exit  for  both  speeds,  in  feet, 

b.  the  drag  of  two  nacelles  for  both  speeds,  in  pounds, 

c.  the  exit  area  of  two  diff users  (one  for  each 
nacelle),  in  square  feet, 

d.  the  weight  of  two  nacelles,  in  pounds, 

e.  the  location  of  the  vertical  center  of  gravity  of 
the  nacelles  below  the  keel  of  the  craft,  in  feet, 
and 

f .  the  location  of  the  longitudinal  center  of  gravity 
of  the  nacelles  forward  of  the  transom,  in  feet. 
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Chapter  7. 

OFF-DESIGN  PERFORMANCE 

Prediction  of  off-design  performance  may  be  carried 
out  by  SUBROUTINE  NACEL  only  after  a  design  has  been  com- 
pleted and  returned  to  the  calling  program.   The  calling 
program  must  reset  the  value  of  a  control  variable  named 
ISTRT  from  1  which  indicates  design,  to  3  indicating 

evaluation,  and  the  value  of  a  control  variable  named 

-9 
DELTA  which  has  a  finite  value  greater  than  1  x  10 

during  design  to  0.  which  prevents  rejection  and  return 

when  cavitation  limits  are  violated. 

The  same  initial  values  are  known  as  when  the  design 
process  started,  and  in  addition,  all  values  established 
during  the  design  cycle  are  now  stored  for  use  during 
evaluation. 

The  free  stream  conditions,  total  pressure,  dynamic 
pressure  and  cavitation  number,  are  calcuated  for  each 
speed  to  be  evaluated. 

The  subroutine  now  proceeds  through  the  same  set  of 
calculations  as  for  design  except  that  the  auxiliary  in- 
let, and  diffuser  sizing  sections  are  omitted.  Cavitation 
checks  are  made  in  the  same  manner,  but  the  existing 
parameters  are  not  allowed  to  be  rejected.   Instead  of 
rejection  and  return  a  variable  named  CAV(I,J)  has  its 
value  changed  from  0.  to  1.  to  indicate  that  cavitation 

is  occuring  at  the  evaluation  speed  indicated  by  I  and  at 
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the  location  indicated  by  J.  J=l  on  the  outside  of  the 
forebody,  2  inside  the  lip,  and  3  in  the  elbow. 

On  completion  of  an  evaluation  cycle  the  subroutine 
returns  the  value  of  the  head  loss  in  feet  and  the  drag  in 
pounds  for  two  nacelles  for  each  of  the  off-design  speeds. 
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Chapter  8. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Results  of  the  subroutine's  use  with  the  optimiz- 
ation program  are  discussed  in  reference  1.   In  addition, 
to  test  it  separately,  the  subroutine  was  called  by  a 
test  program  which  stepped  through  values  of  D^/D     from 
0.4  to  0.9  and  values  of  cruise  Vj/V0  from  0.5  to  1.0  at 
three  cruise  speeds,  40,  45  and  50  knots.  Take-off  speed 
was  assumed  to  be  one  half  cruise  speed.   The  required 
flow  rate  at  cruise  was  arbitrarily  set  at  600  cubic 
feet  per  second  with  the  flow  rate  at  take-off  ten  per 
cent  greater.  The  angle  of  attack  was  set  at  two  degrees 
and  depth  of  submergence  at  five  feet. 

Ram  pressure  recovery  at  the  diffuser  exit  and  cav- 
itation limits  as  calculated  from  the  above  input  are 
displayed  in  Figures  IS   through  21.  The  ram  pressure 
recovery  curves  are  unique  to  the  particular  flow  rate 
combination  assumed  for  this  calculation  and  are  present- 
ed as  samples.  The  cavitation  limits  would  be  the  same 
for  any  flow  rate  at  the  same  angle  of  attack  and  depth 
of  submergence. 

Several  explicit  and  implicit  assumptions  were  made 
in  this  study  to  simplify  it  or  to  try  to  bound  a  real- 
istic value.   Further  investigation  might  prove  useful 
in  any  or  all  of  the  following  areas: 

1.   Very  little  information  on  bell  diffusers  that 
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could  be  directly  applied  in  this  study  was  found. 
Conical  diffuser  performance  was  assumed  to  be  a  lower 
bound  on  that  expected  from  a  bell.   Information  on  bell 
diffusers  similar  to  that  in  reference  11  on  other  shapes 
would  allow  realistic  prediction  of  their  performance. 

2.  The  pressure  coefficient  information  used  from 
reference  3  was  calculated  rather  than  measured.  Although 
the  predicted  values  for  simple  cases  (i.e.,  a  sphere) 
was  compared  with  experimental  results,  no  mention  is 
made  of  comparison  of  calculated  and  experimental  data 
for  a  more  complex  shape  such  as  an  inlet  nacelle. 

3.  The  drag  calculation  was  made  with  the  implicit 
assumption  that  the  drag  of  a  nacelle  with  a  nose  inlet 
and  flow  into  that  inlet  is  very  close  to  the  drag  of 
that  same  body  with  no  nose  opening  or  inflow.  Refer- 
ences 15,  16  and  17  support  this  assumption,  but  only  a 
small  number  of  cases  with  low  D./D  are  examined.  There 
is  some  indication  in  reference  16  that  the  range  of 
V^/V  for  which  this  assumption  holds  decreases  as  D./D 
increases. 

4.  The  effect  of  proximity  to  the  free  surface  on 

drag  was  not  included.   For  overall  system  performance 

comparison,  the  omission  of  this  effect  would  probably 

not  invalidate  the  comparison  but  there  could  be  a  slight 

bias  in  favor  of  the  smaller  nacelle  in  the  case  where 

there  is  a  significant  difference  in  size. 
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In  its  present  form,  SUBROUTINE  NACELLE  requires 
interpolation  in  data  tables  several  times.   Some  of  these 
tables  could  be  approximated  by  functions  to  speed  up  the 
operation  and  reduce  computer  costs.  The  data  for  which 
this  could  be  carried  out  include  Table  12  and  Figure  17 • 
The  data  in  Tables  5  through  16  presents  a  more  difficult 
problem  because  the  relationship  is  not  single  valued  in 
all  cases.   Perhaps  separating  it  into  two  or  more  regions 
with  a  different  function  applicable  in  each  region  offers 
a  solution.  Some  unsuccessful  attempts  were  made  to  ex- 
amine possibilities  for  the  large  data  bank. 

One  other  simplification  which  would  decrease  com- 
puter time  would  be  accomplished  by  elimination  of  the 
calculation  of  diffuser  length.   One  of  the  results  of 
the  data  run  mentioned  above  was  that  the  diffuser  length 
calculated  by  iteration  in  the  subroutine  was  nearly 
always  approximately  ninety  per  cent  of  the  mean  of  the 
maximum  and  minimum  allowed  values.   The  equation 

Ld  =  0.9(Ldmin  +  Ldmax^/2 
could  be  used  to  completely  eliminate  the  iteration. 
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Table  1.   Geometric  Configuration  of  Kuchemann  and  Webber 
Circular  Intakes  Designated  Class  A,  B  and  C 

*  Dimensions  of  all  three  classes  are  defined  by  the 
following  equations  (adapted  from  ref.  7,  p.  79) • 

Do/Di  =  Kx 

Dm/2L  -   K2  -  K3  (Ai/Am)3 

Di/2Li  -   1/K4  (Ki  -  1). 
Table  l.a  Values  of  the  Constants. 


Class 

*1 

K2 

K3 

K4 

A 

1.15 

0.2 

12.5 

1.5 

B 

1.10 

0.2 

12.5 

1.3 

C 

1.03 

0.2 

7.5 

1.0 

(ref.  7) 

Table  l.b  Values  of  the  Primary  Ratios. 


Di/D 


m 


L/Dm 


Class  A  and  B 

Class  C 

0.90 

0.0732 

0.1195 

o.so 

0.1439 

0.231 

0.70 

0.309 

0.452 

0.60 

0.639 

0.910 

0.50 

1.263 

1.575 

0.40 

2.00 

2.175 

0.20 

2.490 

2.494 

*(See  Fig.  3  for  illustration  of  dimension  symbols.) 
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Table  2.  Geometric  Configuration  of  the  Lockheed  Intakes 

Nacelle  *Do/Dm     L/Li     Di/Dm      L/Dm 

A-40-550  0.40      9.0     0.333      2.20 

A-SO-40  0.80      9.0     0.77      0.32 

A-50-280  0.50     9.0     0.426      1.40 

A-75-60  0.75      9.0     0.713      0.45 

A-70-90  0.70     9.0     O.656     O.63 

A-55-280  0.55      9.0     0.433      1.54 

A-60-135  0.60     9.0     0.541     0.81 

A-45-400  0.45      9.0     O.369      1.80 

A-60-200  0.60      9.0     0.541      1.20 

(ref.  3) 

*(See  Fig.  3  for  illustration  of  dimension  symbols.) 
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Table  3»  Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A -40-550 
Internal  Flow: 


Minimum 

Cp 

Vi/Vo 

«=o 

2° 

4° 

6° 

0.70 

0.420 

0.415 

0.400 

0.390 

0.80 

0.260 

0.250 

0.240 

0.220 

0.90 

0.080 

0.050 

0.045 

0.030 

1.05 

-0.260 

-0.255 

-0.285 

-0.550 

1.15 

-0.500 

-0.510 

-0.550 

-1.040 

1.25 

-0.760 

-0.765 

-1.050 

-1.630 

External  Flow: 

0.70 

-0.090 

-0.110 

-0.250 

-0.440 

0.80 

-0.090 

-0.090 

-0.160 

-0.320 

0.90 

-0.090 

-0.090 

-0.120 

-0.230 

1.05 

-0.090 

-0.090 

-0.115 

-0.155 

1.15 

-0.090 

-0.090 

-0.115 

-0.145 

1.25 

-0.090 

-0.090 

-0.115 

-0.140 

(ref.  3) 
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Table  4.   Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-gO-40 
Internal  Flow: 


Minimum  Cp 

Vi/Vo 

cx=0 

2° 

4° 

6° 

0.70 

0.200 

0.050 

O.SO 

-0.100 

-0.235 

0.90 

-0.465 

-0.700 

1.05 

-1.160 

-1.425 

1.15 

-1.7S0 

-2.220 

1.25 

-2.500 

^-2.50 

External  Flow: 

0.70 

-0.435 

-0.430 

-0.540 

-0.60 

0.30 

-0.420 

-0.460 

-0.500 

-O.56O 

0.90 

-0.405 

-0.450 

-0.490 

-0.515 

1.05 

-0.375 

-0.415 

-0.440 

-0.475 

1.15 

-O.365 

-0.400 

-0.425 

-0.460 

1.25 

-0.350 

-0.375 

-0.410 

-0.440 

(ref.3) 
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Table  5.  Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-50-230. 
Internal  Flow: 


Minimum  Cp 

Vi/Vo 

oc=o 

20 

4° 

0.70 

0.3S0 

0.355 

0.320 

0.30 

0.200 

0.160 

0.120 

0.90 

-0.015 

-0.050 

-0.100 

1.05 

-0.370 

-0.420 

-0.435 

1.15 

-0.640 

-0.700 

-0.770 

1.25 

-0.940 

-1.00 

-1.030 

xternal  Flow: 

0.70 

-0.150 

-0.160 

-0.240 

0.30 

-0.150 

-0.155 

-0.200 

0.90 

-0.150 

-0.155 

-0.190 

1.05 

-0.150 

-0.155 

-0.175 

1.15 

-0.150 

-0.155 

-0.175 

1.25 

-0.150 

-0.155 

-0.175 

(ref.  3) 
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Table  6.  Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-75-60 
Internal  Flow: 


6° 


Minimum  Cp 

Vi/Vo 

OCZQ 

2° 

4° 

0.70 

0.250 

0.130 

-0.020 

o.eo 

-0.020 

-0.175 

-0.350 

0.90 

-0.340 

-0.530 

-0.730 

1.05 

-0.930 

-1.230 

-1.600 

1.15 

01.440 

-1.&L5 

-2.370 

1.25 

-2.050 

^-2.50 

^-2.50 

External  Flow: 

0.70 

-0.360 

-0.390 

-0.445 

0.S0 

-0.350 

-0.375 

-0.410 

0.90 

-0.340 

-0.360 

-0.395 

1.05 

-0.320 

-0.345 

-0.370 

1.15 

-0.310 

-0.335 

-0.355 

1.25 

-0.300 

-0.325 

-0.345 

(ref.  3) 
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Table  7.   Variation  of  Minimum  Pressure  Coefficient  with 
velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-70-90. 
Internal  Flow: 


Minimum 

Cp 

Vi/Vo 

£*=0 

2° 

4° 

6° 

0.70 

0.305 

0.215 

0.100 

0.30 

0.070 

-0.040 

-0.130 

0.90 

-0.210 

-0.350 

-0.535 

1.05 

-0.710 

-0.930 

-1.225 

1.15 

-1.130 

-1.420 

-1.920 

1.25 

-1.630 

-2.120 

^-2.50 

External  Flow: 

0.70 

-0.2S0 

-0.310 

-0.375 

-0.455 

0.30 

-0.270 

-0.295 

-0.335 

-0.405 

0.90 

-0.265 

-0.2S5 

-0.315 

-0.365 

1.05 

-0.255 

-0.270 

-0.300 

-0.330 

1.15 

-0.250 

-0.265 

-0.290 

-0.320 

1.25 

-0.240 

-0.260 

-0.230 

-0.310 

(ref.  3) 


-52- 


Table  3.   Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-55-230 
Internal  Flow: 


Minimum 

Cp 

Vi/Vo 

0C=O 

20 

4° 

60 

0.70 

0.415 

0.335 

0.350 

0.30 

0.235 

0.200 

0.160 

0.90 

0.030 

-0.010 

-0.070 

1.05 

-0.330 

-0.375 

-0.440 

1.15 

-0.595 

-0.650 

-0.730 

1.25 

-0.330 

-0.950 

-1.330 

External  Flow: 

0.70 

-0.125 

-0.225 

-0.450 

-0.700 

0.30 

-0.125 

-0.145 

-0.270 

-0.500 

0.90 

-0.125 

-0.135 

-0.135 

-0.320 

1.05 

-0.125 

-0.130 

-0.170 

-0.205 

1.15 

-0.125 

-0.130 

-0.150 

-0.190 

1.25 

-0.125 

-0.130 

-0.150 

-0.170 

(ref.  3) 
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Table  9»   Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-60-135 
Internal  Flow: 


Minimum  Cp 

Vi/Vo 

ex  zo 

20 

4° 

6° 

0.70 

0.300 

0.230 

0.140 

0.040 

0.30 

0.030 

-0.010 

-0.115 

-0.235 

0.90 

-0.175 

-0.230 

-0.405 

-O.56O 

1.05 

-0.620 

-0.755 

-0.910 

-1.135 

1.15 

-0.965 

-1.120 

-1.330 

-1.620 

1.25 

-1.340 

-1.535 

^--1.30 

-2.230 

External  Flow: 

0.70 

-0.255 

-0.230 

-0.315 

-0.330 

o.3o 

-0.255 

-0.275 

-0.300 

-0.350 

0.90 

-0.250 

-0.270 

-0.295 

-0.330 

1.05 

-0.245 

-0.265 

-0.230 

-0.310 

1.15 

-0.240 

-0.265 

-0.230 

-0.300 

1.25 

-0.240 

-0.260 

-0.275 

-0.295 

(ref.  3) 
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Table  10.  Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-45-400 
Internal  Flow: 


6° 


Minimum 

Cp 

Vi/Vo 

<*=o 

2° 

4° 

0.70 

0.400 

0.335 

0.370 

0.30 

0.230 

0.205 

0.135 

0.90 

0.025 

0.000 

-0.020 

1.05 

-0.320 

-0.340 

-0.365 

1.15 

-0.570 

-0.600 

-0.630 

1.25 

-0.350 

-0.330 

-0.915 

External  Flow: 

0.70 

-0.120 

-0.140 

-0.310 

0.30 

-0.120 

-0.130 

-0.210 

0.90 

-0.120 

-0.130 

-0.155 

1.05 

-0.120 

-0.130 

-0.150 

1.15 

-0.120 

-0.130 

-0.150 

1.25 

-0.120 

-0.130 

-0.150 

(ref.  3) 
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Table  11.  Variation  of  Minimum  Pressure  Coefficient  with 
Velocity  Ratio  and  Angle  of  Attack  for  Lockheed  Nacelle 
A-60-200 
Internal  Flow: 


Minimum  Cp 

Vi/Vo 

<x=o 

2° 

4° 

60 

0.70 

0.400 

0.360 

0.300 

0.30 

0.205 

0.150 

0.075 

0.90 

-0.015 

-0.030 

-0.165 

1.05 

-0.400 

-0.435 

-0.620 

1-15 

-0.635 

-0.300 

-1.140 

1.25 

-1.000 

-1.265 

-1.390 

External  Flow: 

0.70 

-0.160 

-0.200 

-0.335 

-0.660 

0.30 

-0.160 

-0.175 

-0.240 

-0.420 

0.90 

-0.160 

-0.160 

-0.210 

-0.230 

1.05 

-0.160 

-0.160 

-0.135 

-0.230 

1.15 

-0.160 

-0.160 

-0.130 

-0.215 

1.25 

-0.160 

-0.160 

-0.130 

-0.210 

(ref.  3) 
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Table  12.  Pressure  Loss  Coefficient  for  an  Elliptical 
Profile  Lip 


Vi/Vo 

AP/qo 

0.6 

0.025 

0.8 

0.033 

1.0 

0.043 

1.2 

0.035 

1.4 

0.130 

1.6 

0.173 

Reynold's  number*^ 

X  10^ 

(Adapted 

from  Fig. 

74,   Report  No.    5,   ref.    3) 
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Figure  4.  L/Dm  vs.  Di/Dm  for  the  Lockheed  intakes  and 
Kuchemann  and  Weber  intakes  Classes  A,  B  and  C. 
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Figure  12.  Minimum  Cp  vs.  L/Dm  in  the  internal  flow  at 
Vi/Vo  =  .80 

-67- 


0 


-.5 
Cp 

-1.0 


-1.5 


1.0       1.5 
L/Dm 


2.0 


Figure  13.  Minimum  Cp  vs.  L/Dm  in  the  internal  flow  at 
Vi/Vo  =  .90 
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Figure  14»  Minimum  Cp  vs.  L/Dm  in  the  internal  flow  at 

Vi/Vo  =1.05 
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Figure  15.  Minimum  Cp  vs.  L/Dm  in  the  internal  flow  at 
Vi/Vo  =1.15 
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Figure  16.  Minimum  Cp  vs.  L/Dm  in  the  internal  flow  at 
Vi/Vo  =  1.25 
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Figure  17.  Ram  Pressure  Recovery  vs.  V^/V  for  the  elliptic- 
al profile  lip  (Blackaby  and  Watson  13E). 
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Figure  18.  Ram  Pressure  Recovery  at  the  diffuser  exit  vs. 
V^/V  for  a  40  knot  cruise  speed 
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Figure  19.  Ram  Pressure  Recovery  at  the  diffuser  exit  vs 
V^/V  for  a  45  knot  cruise  speed 
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Figure  20.   Ram  Pressure  Recovery  at  the  diffuser  exit  vs 
V./V  for  a  50  knot  cruise  speed 
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Figure  21.  Cavitation  limits  as  a  function  of  D./D  vs. 
to  i'  m 

V./V  at  cruise  speeds  of  40,  45  and  50  knots 
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Appendix  A.   CAVITATION  CRITERA 

The  total  pressure  head  at  any  point  in  the  system  can 
be  determined  with  respect  to  the  free  stream  head  by  trac- 
ing along  a  stream  line  from  out  in  the  free  stream  to 
the  point  of  interest. 

H  =  Ho  ~  Hloss 
or  P  =  PQ-  P^oss  i-n  terms  of  pressure. 
Using  Bernoulli's  equation  and  allowing  the  static 
head  term  to  reflect  all  changes  in  elevation  or 
submergence, 

h0  -  V2/2S  --   h  -  V2/2g  -  Hloss 
or  pQ  <-  i/*V2  =  p  ♦  J/»V2  *  Ploss. 
The  cavitation  number,  <r t   can  be  defined  as 

<T=  (Pref.  -  P)/i/9T2ref, 

Where  P .-,  and  V «  are  the  reference  pressure  and 

ref .      rei •  r 

velocity  to  which  the  cavitation  number  is  related  and 
p  is  the  local  static  pressure  at  the  point  of  interest. 

If  the  free  stream  conditions  are  used  for  reference, 
then  p0  and  VQ  are  the  reference  pressure  and  velocity. 
In  this  case 

p0  -  p  ■  i    v2  -  J    v2  ♦  p1oss 
so  <r0  --  (Po  -  pl/i/'v*  =  v2/v2  -  1  ♦  Ploss/i/*v2 

Cavitation  is  assumed  to  occur  when  the  local  static 
pressure  is  equal  to  the  vapor  pressure  of  water.  The 
incipient  free  stream  cavitation  number  is  defined  by 

^oi  =  (P0  -  Pt>/*/"S  =  <Vcrit/  *l]-1   +  Ploss/^Vo- 
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This  definition  allows  the  prediction  of  cavitation  using 
the  free  stream  incipient  cavitation  number  if  the  losses 
are  known. 

Consider  the  situation  in  Figure  (3).  The  diffuser 
entry  is  designated  station  1,  and  the  diffuser  exit 
station  2,  to  remain  consistant  with  the  rest  of  the  text. 
The  turning  vane  is  not  intended  to  be  a  particular  type 
or  design,  but  is  expected  to  be  well  designed  and  sized 
for  the  elbow  and  diffuser  with  which  it  is  used.   The 
sketch  is  conceptual,  only.   The  subscript  tv  will  designate 
conditions  on  the  turning  vane. 

To  determine  the  cavitation  number  on  the  turning  vane, 
conditions  at  station  2  can  be  used  and,  in  turn,  related 
to  the  free  stream.  As  before, 

Ptv  =  P2  "  Ploss 

2 — ►tv 

If  the  turning  vane  starts  near  station  2,  then  any 

loss  between  there  and  the  point  of  lowest  pressure  can 

be  neglected  as  it  will  be  very  close  to  the  leading  edge. 

Ptv  =  P2  =  Ptv  +  iKv  =  P2  +  ^V2 
The  cavitation  number  on  the  turning  vane  is, 

<Ttv=  <P2  -  V*fvr 

When  ptv  =  py,  cavitation  starts,  so  the  incipient 
cavitation  number  is 

r2 


tftvi  =  (P2  -  PvVi/*v| 


^~tvi   *s   a   characteristic  number  determined  from 
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P2  =  *~tvi  */°i  +  PV 


experiments  and  is  therefore  known  (ref.  28).   The  con- 
ditions at  station  2  which  will  cause  cavitation  on  the 
turning  vane  can  be  determined.   Solving  the  equation 
above  for  p2, 

P2  =  *"tvi 

but  P2  can  be  determined  in  terms  of  the  free  stream 

conditions,  which  are  known. 

p     =  p     +,   1/7  v2   -   h/>Vo   -   P, 
v2       *o  '      o        */     2  loss 

or     tfl    •    i/'V2   +  P     =  Prt   +  \PV2   -   |/>V2   -P., 

wtvi  A/     2       rv   ^o   ^a  0   *r    2         loss 

P,    in  this  expression  is  the  sum  of  all  losses  from  the 
free  stream  to  station  2.  V  is  the  only  variable  re- 
maining which  indicates  the  conditions  at  station  2,  so 

2 

solving  for  V  , 

v*  (1  +  <r  ,)  =  f(p  -  p  )/i/*v2  +  l  -  P,   /i/>v2]  v2 

2      "tvi     I  *o   rv  /*/     0        loss   '   oJ  o 

The  first  term  in  the  brackets  is  the  free  stream 

incipient  cavitation  number.   The  value  of  V_  which  gives 

rise  to  incipient  cavitation  on  the  turning  vane  is 

V0  =    f  n-.    +   1  -   P.        /MV2   1*     V  /    (1  4     rr     .)*. 
2        L^oi  loss' ^  /     o  J  o'  utv± 
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Appendix  B.   GEOMETRIC  DERIVATIONS 

There  are  three  variable  lengths  in  the  nacelle  which 
are  not  a  function  of  the  forebody  geometry  and  which  must 
be  calculated  in  each  design.   They  are  the  diffuser  lengta , 
L  ,  the  additional  length  due  to  the  auxiliary  inlet,  Ly, 
and  a  length  after  the  diffuser  exit  containing  part  of 
the  elbow  which  directs  the  flow  upward  into  the  strut,  the 
strut  inlet  length. 

1.   The  Strut  Inlet  Length 

Assume  that  the  thickness  to  chord  ratio  of  the  strut 
all  along  its  length  from  the  nacelle  to  the  free  surface 
is  determined  by  the  criterion 

tm/c0  =  [(1  +  Ol)i-   l]/l.l5       (ref.  1,4) 
where  G\   is  the  local  free  stream  incipient  cavitation 
number.   Also  assume  an  elliptical  cross  section  for  the 
strut  and  a  rectangular  cross  section  for  the  internal 
ducting  (ref.  1). 

The  respective  areas  are 

Astrut  =  {7TA}  tm  co 
Aduct  a  2  tm  c0  (ref.  1) 

Aduct  mus<t  be  tne  same  as  the  exit  area  of  the  diffuser, 

A2- 

To  allow  for  structure  and  fairing,  assume  that  the 
desired  additional  nacelle  length  is  close  to  the  value  of 
the  strut  chord,  c0. 

For  0~±   =  0.3,  tm  =  0.12  cQ,  2A2  =  0.12  cQ2,  so  that 
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c0  =  3.545  D2   . 

2.  Length  Due  to  the  Auxiliary  Inlet 

If  the  auxiliary  inlet  has  the  configuration  shown 
in  Figure  3,  the  inlet  area  is  a  section  of  a  right  circ- 
ular cone  with  a  minimum  diameter  D.,  and  slant  height  X. 

The  angle  of  inclination  from  the  centerline  of  that 
cone  section,  shown  in  the  figure  as  0X,  is  assumed  to  be 
formed  between  the  centerline  and  a  line  passing  through  the 
nose  leading  edge  to  the  point  of  tangency  of  the  forebody 
at  maximum  diameter,  Dm. 

9X  =  tan-l  [(Dm  -  Di)/2L] 

The  mean  diameter  of  the  conical  section  is  then 
D  =  D±   +  X  sin  6X 
and  its  area  is 

Acone  =TTDX  • 
The  area  of  the  inlet  will  actually  be  some  percentage 

of  the  area  above  because  some  allowance  must  be  made  for 

structure. 

AAux  =  K  TTD  X 

K  =  O.S  is  used  in  the  subroutine. 

Finally,  solving  for  X  in  the  above  expression, 

X  =    (0.5/sin  0X)    [(D-i2   +   1.273   AAux  sin  6X/K)^-   DjJ    , 

and  Lx  =  X/cos  6X     . 

3.  Diffuser  Length 

The  length  of  the  diffuser,  L,,  becomes  a  factor  in 

the  determination  of  nacelle  length  when  the  sum  of  the 
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LN  =  LLip  +  Lx  +  Ld  +  3.545  D2  -  5-5 


internal  lengths  exceeds  the  minimum  length,  5»5Dm»   The 
excess  length  is  chargeable  to  the  diffuser  because  all  the 
other  internal  length  requirements  are  fixed  by  other  con- 
siderations.  The  difference  in  length  is 

Dr 
The  change  in  drag  is  given  by 

ADrag  =  CD  |/>VQ2  AS, 
because  CD  is  essentially  a  constant  for  small  changes  in 
LN. 

The  change  in  length  can  be  considered  to  be  all  in  a 
cylinderical  section  of  diameter  D  • 

AS  =  77 Dm  ALN 
The  increase  in  power  due  to  the  increase  in  drag  is 

APower-]_  =  CD  |/7V0377Dm  A  LN   . 
The  power  to  overcome  the  diffuser  losses  is  found  by 

APower2  =  (CKt  +  fLd/dm)J/>Vi2Q   . 
The  total  increase  in  power  required  to  operate  the 
diffuser  is  then 

Al  Power  =  APower-^  +  A  Power2 
In  the  subroutine,  the  minimum  value  of  A  Power  is 
determined  by  iteration  on  Lj. 


-82- 


Appendix  C.   SUBROUTINE  LISTING 

SUBROUTINE  NACEL  and  its  flow  diagram  listed  in  the 
following  pages  are  in  the  form  for  which  the  rationale  is 
described  in  chapters  6  and  7.  Similar  listings  for 
FUNCTION  TABLE,  which  was  not  programmed  by  the  author  but 
was  furnished  for  common  ues  with  reference  1,  follow 
NACEL.  FUNCTION  CFS  is  also  shown  for  completeness  but  no 
flow  diagram  was  generated  for  it  because  of  its  simplicity. 

The  FORTRAN  names  below,  used  for  various  variables  and 
constants  in  the  subroutine,  are  described  as  an  aid  in 
understanding  the  logic  and  specific  calculations.   The 
symbols  are  listed  in  approximate  order  of  first  occurance 
beginning  with  the  first  executable  statement. 

ZK  -      A  constant  decimal  less  than  one  which  in- 
dicates how  much  of  the  annulus  occupied  by 
the  auxiliary  inlet  is  available  for  opening. 

SPO  -     Static  pressure  in  the  free  stream  in  p.s.f. 

HS  -      Depth  of  submergence  in  feet. 

HA  -      Pressure  head  of  the  atmosphere  in  feet  of 
water. 

RHOW  -    Density  of  water  in  slugs  per  cubic  foot. 

G  -      Acceleration  of  gravity,  32.174  feet/sec? 

PVP  -     Vapor  pressure  of  water  in  p.s.f. 

PV  -      Vapor  pressure  of  water  in  feet  of  water. 

VELR(I)  -  Inlet  velocity  ratio  at  speed  I.  1=1 

indicates  cruise.  1=2  indicates  take-off. 
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I  greater  than  2  indicates  a  performance 
evaluation  speed. 

VI(I)  -   Inlet  velocity  in  feet/sec. 

V0(I)  -   Free  stream  velocity  in  feet/sec. 

SIGTV  -  Incipient  cavitation  number  on  the  turning 
vanes  referenced  to  diffuser  exit  pressure 
and  velocity. 

JNUMB  -   A  counter. 

ISTRT  -  A  control  variable  indicating  whether  the 
current  calculation  is  part  of  the  design 
process  (ISTRT  =1),  or  part  of  the  eval- 
uation process  (ISTRT  =3). 

NUMB  -    A  control  variable  indicating  the  total 
number  of  speeds  to  be  considered. 

QO(I)  -  Free  stream  dynamic  pressure  at  speed  I  in 
p.s.f . 

SIGI(I)  -  Free  stream  incipient  cavitation  number. 

PTO(I)  -  Free  stream  stagnation  pressure  in  p.s.f. 

TRIM(I)  -  Angle  of  attack  in  degrees. 

CPEX  -    External  minimum  pressure  coefficient. 

CDUMX  -   A  dummy  array  generated  during  interpolation. 

CDUMY  -   A  dummy  array  generated  during  interpolation. 

CPEXT(J,K,I)  -  Tabulated  minimum  pressure  coefficients 

in  the  external  flow  at  the  Jth  velocity 

ratio,  Kth  length  to  diameter  ratio,  and  Ith 

angle  of  attack. 
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CPINT(J,K,I)  -  Tabulated  minimum  pressure  coefficients 

in  the  internal  flow. 
VRT  -     Tabulated  array  of  velocity  ratios. 
IL,  KL,  ML,  NL  -  Control  variables  indicating  the 

order  of  interpolation  desired. 
XDT  -     Tabulated  array  of  angle  of  attack. 
XDTT(I)  -  A  dummy  array  generated  during  interpolation. 
DIDMT  -   Tabulated  array  of  Di/Dm. 
DIDMX  -   The  maximum  value  of  B±/Dm   for  a  Vi/V0  as 

determined  from  interpolation  in  the  data 

table. 
XD  -      A  working  variable  for  length  to  diameter 

ratio. 
WGTS(1,1)  -  Wet  weight  of  the  nacelle  in  pounds. 
DIDM  -    Current  value  of  Di/Dm. 
QI  -      Flow  rate  at  cruise  through  one  inlet  in 

c.f.s. 

Nose  inlet  area  in  feet  squared. 


AI  - 
DI  - 
DM  - 
ELEXT  - 


Current  value  of  D^   in  feet. 
Current  value  of  Dm  in  feet. 


Length  of  the  forebody  in  feet. 
VRTT(I)  -  A  dummy  velocity  ratio  array  generated  during 

interpolation. 
VRTEX(I)  -  A  dummy  velocity  ratio  array  generated  during 

interpolation. 

VRMAX(I)  -  Maximum  allowed  velocity  ratio  for  a  given 
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D-/Dm  as  determined  from  the  tables. 

QIN  -     Working  variable  for  flow  rate  through  the 
nose  inlet  at  take-off  in  c.f.s. 

QC  -      Required  flow  rate  at  take-off  in  c.f.s. 

QAUX  -    Flow  rate  through  the  auxiliary  inlet  in 
c.f.s. 

KDEX  -    A  counter. 

VI2  -     Working  variable  for  VI(2). 

PRLT  -    Tabulated  array  of  pressure  recovery  coeffic- 
ients for  the  lip. 

PRLT2  -   Lip  pressure  recovery  coefficient  at  take- 
off. 

PTI  -     Stagnation  pressure  in  the  internal  flow  in 
p.s.f • 

SPI  -     Static  pressure  in  the  internal  flow  in  p.s.f. 

AIAUX  -   Area  of  the  auxiliary  inlet  in  feet  squared. 

VTAUX  -   Auxiliary  inlet  velocity  ratio. 

PTAUX  -   Stagnation  pressure  inside  the  auxiliary 
inlet  in  p.s.f. 

PRAUX  -   Pressure  recovery  coefficient  of  the  aux- 
iliary inlet. 

DYP  -     Dynamic  pressure  inside  the  auxiliary  inlet 
in  p.s.f. 

PC  -      Stagnation  pressure  in  the  combined  flows 

from  the  auxiliary  and  nose  inlets  in  p.s.f. 

DLIP(I)  -  Inlet  losses  at  speed  I  in  p.s.f. 
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QDIF  -    Flow  rate  through  the  diffuser  in  c.f.s. 

ELENT  -   Interior  length  of  the  lip  in  feet. 

PHI  -     The  angle  labeled  6X  in  Figure  3. 

ELAUX  -   Interior  length  due  to  the  auxiliary  inlet 

in  feet. 
ELMAX  -   Maximum  permissible  length  of  the  diffuser 

in  feet. 
ELMIN  -   Minimum  permissible  length  of  the  diffuser 

in  feet. 
QDIF  -    Flow  rate  through  the  diffuser  in  c.f.s. 
EL  -      Current  value  of  the  diffuser  length  in  feet. 
DEL  -     Change  in  diffuser  length  in  feet. 
ELN  -     Length  of  the  nacelle  in  feet. 
ELD,  ELL  -  Working  variables  for  ELN. 
DDM  -     Mean  diameter  of  the  diffuser  in  feet. 
XKT  -     Sudden  expansion  loss  coefficient,  Kt. 
REL  -     Reynolds  number  based  on  length. 
RED  -     Reynolds  number  based  on  diameter. 
DL  -      Ratio,  Dm/L^. 
CDRG  -    Drag  coefficient. 
ANGL  -    Deffuser  half  angle  of  expansion,  labeled 

6  in  Figure  3 • 
POW  -     Current  value  of  the  power  charged  to  the 

diffuser  in  ft. lbs. /sec. 
POWI  -    Previous  value  of  POW. 

ELDIF-    Final  value  of  diffuser  length  in  feet. 
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REND  -    Same  as  RED. 

RENL  -    Same  as  REL. 

DDIF  -    Diffuser  loss  in  p.s.f. 

PLOSS  -   Sum  of  diffuser,  pipe  and  lip  losses  in  p.s.f. 

VAOUT  -   Average  diffuser  exit  velocity  in  ft. /sec. 

VCRIT  -   Critical  diffuser  exit  velocity  at  which 

cavitation  occurs  on  the  turning  vanes  in 

ft. /sec. 
VMAX  -    Maximum  diffuser  exit  velocity  in  ft. /sec. 
DELH(I,1)  -  PLOSS  converted  to  head  in  feet  of  water 

at  speed  I. 
PRC(I)  -  Ram  pressure  recovery  coefficient  at  the  dif- 
fuser exit. 
CD(I)  -   Drag  coefficient. 

AEXN  -    Wetted  area  of  the  nacelle  in  feet  squared. 
CPOD(I)  -  Drag  of  two  nacelles  in  pounds. 
AREA(l)  -  Diffuser  exit  area  of  two  diffusers  in  feet 

squared. 
CGS(1,1)  -  Distance  of  the  vertical  center  of  gravity 

of  the  nacelles  from  the  keel  in  feet. 
CGS(2,1)  -  Distance  of  the  longitudinal  center  of 

gravity  of  the  nacelles  from  the  transom  in 

feet. 
CGS(3,1)  -  Distance  of  the  VCG  of  the  contained  water 

from  the  keel  in  feet. 

CGS(4,D  -  Distance  of  the  LCG  of  the  contained  water 
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from  the  transom  in  feet. 
RHOD  -    Density  of  steel  in  slugs  per  cubic  foot 
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cha^t   Tin  r   -  supr-iut  in 


VITOFLU     ft- 


IK  IS    THE    OEC I" 

PART  IF  THC  ANNUL 
OCCUPIED  BY  THf 
AUKILlAkY  INLET  T 
IS  ACTUALLY  OPFNI 
THE  BEMAJNnER  IS 
STRUCTURE. 


1 

CI 

1                    IK     •     .8 

1          spo  .    ihs  , 

I                H»|»»MCW*G 

1      pvp  *   »v«RHrb» 

1                 VFl«(l 1    • 

1             VII1I/VGI1I 

I 
I 

SIGTV   IS  THE 
INCIPIENT  CAVITATION 
HO.  ON  THE  ELBOW 
TURNING 

VANES  REFERENCED  TP 
OIFFIISER  EXIT 
PRESSURF  ANO 
VELOCITY. 

I 

I 


1             SICTV    - 

0.* 

1                 JNIJM8 

-    2 

NUMB       IS    AN    INOEI 

INDICATING     THE     TOTAL 

NO.    OF     SPFFDS    TO    flf 

EXAMINED. 

NLH8*2     IF     ONLY    CRUISE 

ANO    TAKE-OFF    ARE 

SPECIFIFO. 

NUMB*  3    0«     "     -'■      IF    ONE 

0«    MORE    OFF-CES  IGN 

SPEEDS     ARE     SPECIFIEO. 

I 

I 

I 
1  STUf        IS     AN     INDEX 
INDICATING    yHFTHFR 
THE    CURRENT    OPERATION 


PART    OF    THE 

DESIGN 

PROCESS, 

IN 

MH 

CM    ITS 

VALUE     IS 

1. 

OR 

PART 

OF 

THE    EVALUATION 

PROCESS. 

IN 

MH 

CH 

CASE     ITS 

VALUE 
1 

IS    3     . 

ISTRT     .EC.    3 


ITPUf 

I 
I 
I 


I       NOTF       05 


8FGIN   or   LOOP 
10    I    «     ISTRT, 

JNUKB 


-1 

. .s* — >• 

1                       04 

1                  0" (  1  1     '                  1 

1             ,S,HH  1N"V   Mil             | 

I            •vim              i 

I 

I 

Gl III    IS   TMt 
riPIENT    CAVI TATION 
.     "FFCRENCED    TO 
Ec     STREAK 

iditions. 


TRI«(I I    .GT. 


TPI«M 11-3. 
I 


I 
.    I 
INTERPOLATE     IN    TMF 
DATA    TABLE    TO    FINO 
THE     INLET    KITH    THE 
0ES1RFD 
PRESSURE    COEFFICIENT. 

I 

I 


1 

NOTF 

12 

. 

BEG! 

N    00 

LOOP 

. 

* 

610 

,     2 

* 

01 

,?0  — 

1 

->l 

1 

NOTE 

1) 

-01.  !<,.-->• 

I  IS 

I  CDUHIIJI    ■  I 

I         CPE'TU.K.II  I 

i 

I 

608  •         16 

NO       •  • 

--•        ENO   OF    on        • 
•         LOOP?         • 


1 

17 

CDUMVIK 

TABlflvRT, 

VFLBI1  1 

>     ■                 1 

cdumk,     i 

,NLt              1 

XOTTIII     - 
ABIE.CDUMV.XPT, 

f  a  i, m  I 


I     TABLFI XDT.OIOMT. 


I 

I 
IF     THE     TRIAL     NACELLE 
HAS    LESS    FRONTAL     ARFA 
THAN    THE     MINIMUM 
REQUIRED     TO    AVOID 
CAVITATION,     REJFCT 
THE    TRIAL     VALUE.        IF 
NON- 

CAVITATING,     CALCULATE 
INLET    DIMENSIONS. 


I 


I     WGTSII .11     »     l.E 13     I 


NOTE      23 

•             OIDNK 

.IT.               • 

•    010N*ll. 

»    OELTAI    • 

•    .ANO.    DELTA    .GT.       • 

»                l.E 

-    9                 • 

-101- 


CHART     till  i 


»  tl'TINt        *Alcl 


(H»l     Stt    - 


11  .Z«—  ->l 


I  0!     -     .S.3ILI 

I 

I       ai  *    ji/vini 

I 

I  11   = 

I     SORT  (All  «l.l2t-M 

I 

I  P«    '    "II  )I3M 


flfTFOMINE  MAI.  CI  )M 
PATF  A»  T4K;-0fC  ANO 
=  JUI»L0 


Fir 


AIE. 


t'l 


•  ikiiiart  inlft  must 
p.f  si7c0  tp  accept 
any   excess 

PE0U1SCD  Ft_TK. 


I    FLEXT  «  0N«I0 

I 

I  CPIN  «  -  SIGH ?l 


I 

INTEFPOl  ATF      IN     TMC 

TATA     TABLE     TO 

DETERMINE    THE    «»IMIt 

VELOCITY 

PATJCS    AT    CRUI St     ANT 

TAKE-TFF . 


I 

I 

I       NDTF        13 


I 
I 

ASSUMING    THE    aux. 
1NLFTS    ALIO*    FLOW    TO 
ENTER    BEFORE     THE 
OIEF'JSFR  f 

CALCULATE  LOSSES  ANO 
TOTAL  PRESSURE  OF  THE 
COMBINED    FLOW. 

I 

I 


1 

IT 

1          DC 

l.S"UI.-'l            1 

1       QAUI 

"  oc 

-  aiN    i 

I        N3TF       1* 


I  COUMfUl     * 

I  CP)hT(»,j,|l 


VRTTIII    ■ 

r*fti.6<CDU*«Y,  V«T, 
CPIN, U   t 

vptfiiii    » 

TABt  E(CDUMT.V->T, 
CPEX.lLl 


1               OAUI    ■    0. 

1 

1 

I     oi n  ■  ac  -  u«u> 

1               K.DEX    ■    0 

Ires 

I 

I 
I 
I 

I 


la 


I  COUNVUI     « 

I     TABLE) XOT.CDUNX, 

I  XO.Kll 

a. . 

I 
1 

TP<J  •         09 

•  • 
•          FNP   OF    03 

•  LO'IP?  • 


IVES 

I 
I 
I 
I 
I 


VRMAXllI     >  | 

TABLFIALFAT.  I 

V°T(   ■  ,T»  |M|  1    I.HL I  I 

I 

VRMAXI2I      ■  | 

TABLE IALFAT.VRTT.  I 

TRINm.Mll  I 


I 

CHECK    FOR    LIP 
CAVITATION    AT     CRUISE. 
IF    tuiunv., 
REJECT. 

I 

I 

I   NOTE   13 

•  VELBCIl  .ST.     • 

•  VRMAXllI  .ANO.    • 

•  OFLTA  .GT.  1,6  -   • 

•  Q  • 

I 

I 
I 

•     1* 

•  •  FALSE 
•    SEE  NOTE     • 

•  ABOVE    " 


I 

I 
CALCULATE     STATIC 
PRESSURE     IMMEDIATELY 
AFT.    OF    THE    LIP. 

I 

I 


I  VIZ    -    QIN/At 

I  V«    •    VI2/V0IJI 

I 

I  PRU    " 

I  TABLFIVRT.PRLT. 

I  VR.JLI 

I 

I  PT|    • 

I  PRL2"80I?I    ♦    SPO 


1 

SPI    ■ 

PTI     - 

1    0. 

5«RH0W 

»VI2«VI2 

CALCULATE 

1 
1 
1 

COMBINEO 

FLOW 

PRESSURES    ANO 

AU«. 

INLET 

AREA. 
1 
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CHART    TITLf     -    Svl«»-  IUTINI       NiCil 


I 

»HU«    »  C. 

YIAUX    .  J. 

•TAIIX    »  C. 


I  PRAlJX    a      ). 

I 

I  PTR'JX     ■ 

I     PRAUX«00I2I 

I 

I    DrP    -    PTAUX 

I 


SP} 
S»l 


I 


t  !»■!« 


1     S0BTI2.«irP/fiHC*l 


I 
I 

IMF  TOTAL  PRESSURE  1 
THE  COKRINEO  FlU  15 
CUCdllTFl    As    TMF 

hass  weighted  avr4g 
of   the   combining 

FIOWS. 


I  vim  - 

I  SQPTI'.pIPC    - 

I  SPI  I/RHOHI 

I 

I  01IPI2I    »    I.    - 

I  ipc  -  srai/gim 


1 
1 
1 

17 

1               OOIF    -    3C 

1 

I 
I 

calculate  the 
internal   lengths  to 
the  oiffusfr  entry. 

I 

I 

I 

» 

I     ELENT    =    ELE»T/<>. 

I 

I  PHI     » 

I  ATANI.5*ION    - 

I  OII/EIFIT1 

I 

I         PHI     -     SINIPHI I 


.5«ISORTIOI»A2    ♦ 

1.37324«»IAUX» 
PHS'ZX)    -    II l/PHS 


SUE    THE    0 


I        m  ■  ».w»"         i 
I  I 

I  Dl     >    01  I 


02    .IT.    01 


02-01  I 


FLNIN    ■ 

2.936075M0?    - 

Oil 


OECIOE     WHICH 
CONDITION    GOVERNS    THE 
0IFF1JSER. 


01     .01.    OC 


/       ii       / 

;  • 

1       NOTE 
C 3NTIVUF 

1* 

IT 

i 
i 

11    •    1 
OOIF    ■    01 

1 
1 
1 

1 

1*             1       NOTE 

IS 

• 

CONTINUE 

' 

I  El     »    Fl     •    Del  I 

I  I 

I       «0E»    ■    OEX    *     1       I 


•    XDEX    .GT.    13 


EIO    -    ELENT    . 
El    ♦    3.«>*<.'J1»02 
ElAUX 

Ell     ■    5.5«0H 

EIN    •    EIO 


ITKUE 

I 
I 
I 
I 
I 
I 


I       / 
/    ♦.01 
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CM4H     T|TLr      -     SUPUHJUHf        NACtl 


iUinn.n    CHAkT     SFT     - 


flfac   .Lt.   r . 


Ei^If        IS     THE 

niFfUSfP     LENGTH 
0P3UIP  ING    TMg    I  FAST 
TOT  At     PPWEP 
FPU     THE     OfcS.Rfn 
DIFFUSION     PATID. 


I  TKUF 

I 


flPIF    =    El 


CALCULATE     THE     LIP 
LOSSES    FOB     EACH 
SIT'lATirW*. 


VOI  I  I  l-FLN/GNU 


1     piipcii 

1        TAILFtV 
1             VEiPI 

RT 
11 

1.     - 
.PRLT, 
.  JL  1 

01     =    OM/FLN 


COOS    » 

cfsipeiihi.   < 

1 . S»OL" • I 3/Z I 


FL  .IT.  3.)0l 


I  «NOl  "  »'••<  111  '  -  I 
I  011/12. »El I >  I 
I       •!!.  2958       I 
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I 


I      cniF  • 

I  3.19E-3«»NSL« 

I       »NGL  » 

I  ■.»t?E-t««MSL 

I 
I 
I 

P)4  « 

»HOW«VOI  I  I  I 
•  •3«£LF«C     ♦ 
(OlF.xKf     • 
FRICTIOEDI 
•EL/OOHI 
•.S«»Hng«VI<  I  II 
•VII II !• jrif 

(IF  I.    •    .t«IFLMAX 


*OEI    .13.     1       • 


I 


I 

I  10 

I         OEl    •    IPOUI    -  I 

I  PriMI*Fl/POWt  I 

I 
>l 

m        I  n 

I  powi    •  POX  | 

I 
I 
I 
•      l? 


NUMB    .IT.     5 


IF»ISE 
I 


OLIPIJI     ■    1. 

TiOLEIVBT.PftL 
VFL9I Jl . JLI 


•         ENO   OF    00 

•         IO0PT         ■ 

•         • 

• 

ITES 

1 
1 
1 
1 

11!           1 

19 

1               JNUMB    ■    2 

1 

1 

— • 

ISTPT     .£0. 


IFilSF 

I 
I 
I 
I 
I 
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HUTOFL-h    f  H40T     S£  T 


CH40T     tl'lf     -     ^"<'   HITINf        NAC^l 


1 

JNUMK     a     Nil*! 

1 

L* 

1 
•  JO >  1 

1       NCTE 

C2 

. 

BEGIN    OC    LPnP 

. 

• 

17     1    »     ISTRT, 

» 

" 

JNUH1 

• 

05 

1 

.22 >l 

1 

C3 

1 

C«V( 1,31    »   r. 

1 

I 

I 

calculate  the 

OtFFUSER     ANO    P|PE 
LOSSES    FOR    EACH 
SITUATIC'N    ANO 
•00    TP    TNF     IIP 
LOSSES. 

I 

I 

I 


no  if  ■ 

ICPIFPXKT    ♦ 

FRICTIRFNDt 

•EiniF/COMP 

•  0.-."PMOM«VTU  I 


PLOSS  • 

olipii  i»aom  • 

ooif  ♦ 

frict1renp' 

•elaux/oi'.s* 

«HPW»VI  t  1  I  •  VI  I  1  I 


I 


PLOSS    - 
Ol|P<2l«OOI2l 
00  IF 


VAOUT    ' 

01  II 

•0.636*2/tO2«D?l 


I 

I 
OETERNINE     THE 

emeu    10CU 

VEIOCITY    AT     THE 
MFFUSE"    EXIT    AT 
yHlfM    CAVITATION    ON 
THF    TURNING    VANES 
OCCURS. 

I 

I 


I 


10 


I  VCKIT    •  I 

I     SQRTISOUARI-VOII  I     ) 
I     /SOPTH.    •    SICTVI    I 

I 
I 
I 

ESTIMATE  THE  MAXIMUM 
L1CAL  VELOCITY  AT  THE 
OIFFHSEK  €«IT. 
I 


11 


I     VMAX    »    l.SO*VAJUT     | 


I 
I 

IF    CAVITATION    OCCURS, 
REJECT    ON    OESIGN, 
1*01 CATE    CN 
FVALUATIPN. 

I 

I 


•   *    •    • 

•    • 

•  •  •  •  * 

•       VH»X 

.01 

VC«|T       • 

•     .AND. 

rsiar   .eo.     • 

•           1     . 

iNn. 

DELTA          • 

•          .GT 

,  i.e  -  *       • 

>l 

AT    THIS    POINT    THE 
OIFFOSFP    HAS    8EFN 
SI7E0    TO    AVOIO 
CAVITATION    AT 
BOTH    TAKE/OFF    ANO 
CPUISF.     INTFRNA1     F 
LOSSFS    ARE    DETEFM| 

I 

I 

I 


The  FOLLOWING  CARP  is 
USED  ONir  WITH  A  TEST 
PROGRAM. 


I      PRCII I  "       I 
I  (PTOIII  -  PLOSS  -  I 

I         spai/oom  I 


com  ■ 

CFSIRfNl  tall. 

1.;«<om/eln> 

•  ••3/21    • 

?.P|DM/ELNI«»3 


I 

I 
CALCULATE    WETTEO 
SURFACE     ANO    ORAG. 

I 

I 

I 


.  S0RTI1.  • 
•II2.»ELEXTI 
/0«I««2I 


•f XN    ■    AEXN    *           1 

P|»OM*CELN    -            | 

?.*FifiTj              1 

CP^i !  i   '            l 
?.«ani!i              1 

*AEX*-rD(I»              1 

1               AREAtll     > 

1          PI«02*02»O.S 

1       CCSIlill    ■    HS 

.       | 

1               HE     -   HCL 

1     CGSI2.1I    -     US 

.       | 

1               .5«IELN    - 

1             3.5*«1»D21 

1 
1 
1 

• 

THATAI1I  .GE. 


IFALSE 

I 
I 
I 


THATAIl  1/57. 29S8 

CGSI2.ll  ' 

CGSI2.ll  • 

CGSU.ll/TANIPHII 


>l 

IS 

1 

2 

COS  13 

•  11 

- 

0. 

CGSK. 

.1  1 

- 

0. 

1 

1 

1                       2 

WCTSIl.ll     =■ 

RHOO    -    RHOU)     * 

1S.0T-*AREM U 

•CEtENT     ♦    ELAVJX 

FLI 

I 
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Appendix  D.   ALTERNATE  DESIGN  RATIONALE  USING  ONLY  VELOCITY 
RATIO  AS  A  SYSTEM  PARAMETER 

Another .approach  to  the  design  of  the  inlet  nacelle 
can  be  made  using  the  maximum  allowable  D^/D  at  a  given 
V^/V  in  stead  of  allowing  it  to  vary  as  a  parameter.   The 
design  proceeds  as  previously  described  except  that  the 
diameter  ratio  determined  from  interpolation  in  the  data 
table  is  used  to  calculate  the  nacelle  geometry.   If,  after 
the  diffuser  is  sized,  the  exit  velocity  is  too  high,  then 
instead  of  rejecting  that  set  of  parameters,  a  new  value 
of  Dm  based  on  the  maximum  allowed  exit  velocity  is  cal- 
culated and  the  subroutine  goes  back  to  the  beginning  of 
the  design  phase  where  a  new  design  begins  with  this  value* 

Performance  evaluation  proceeds  exactly  as  before. 

The  following  pages  contain  a  listing  of  the  subroutine 
in  this  form  with  a  flow  diagram. 
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INDICATING 
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TJTAL 

NO.    OF    SPEE 

3S 
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EXAKINFn. 

NU»R=?    IF    " 

JLT 

CRUISE 

•NT    'AKE-PF 

■    APE 

SPECIFIED. 

NIJwB-3    OR    M 

:re 

IF    ONE 

00    MORE    OFF 

-CE 

IGN 

SPEEDS    ARE 
1 

5PECIFIEO. 

1 
IST«»       IS    A 

1    1 

IOEX 

INDICATING    WHETHER 

THF    CURpENT 

*P"  ATI     N 

IS 

PART    OF     THE 

CE" 

ION 

PPX'SS.     IN 

WHICH     ITS 

VALUE    IS     1, 

OR 

PART 

OF 

THF     EVALUATION 

PROCESS.     IN 
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CASE     ITS    VALUE 

IS    3    . 

ISTRT    .[0.    5 
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1 
>| 
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05 
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BEGIN    OC    LP^P 
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JNU"B 
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0" 

*  1   = 

.S»PH 

jy*v 

III 

*v 

■t  !l 

1 
I 

SICK  I  I    IS    THF 
INCIPIENT    CAVITATION 
NO.    REFERENCED    To 
FREE    STREAM 
CONDITIONS. 


I     SIC.  Mil     =    (5P0 
I  PVPI/'JOIII 


TRIX1I 
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SIGII1)        I 


1 
I 

INTERPOLATE    IN   THE 

DATA    TABLE    To    FIND 

THE     IALET    WITH    THE 

OESIPFO 

PRESSURE    COEFFICIENT. 
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I       NOTE       H 


01.  20 >l 


I       NOTE       13 


1* 

1          coiimx 
1        rocxil 

il    » 

1 
1 

I  X0TT1II     • 

I     TA8LFIC0U«V,X0T, 
I  CEX,«LI 


010"    = 

TARLEIXOT.DIOMT, 

xn.FL I 


I 
I 

IF    THF    TRIAL    NACELLE 
HAS    LESS    FRONTAL     AREA 
THAN     THF     MINIMUM 
REQUIRED    TO    AVOID 
CAVITATION.     REJECT 
THE    T=IAL    VALUE.       IF 
NON- 

CAVITATING.    CALCULATE 
INLET    DIMENSIONS. 
I 
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OECX    FOB    LIP 
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ISE. 
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Pi 

JECT. 
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NHTf 
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■     ■ 
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G». 
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V»»AXUI     . 

AND. 
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>« 

DET£P*IN?     HAX.     =  L"W 
PATC     AT     TAK=-fFF     AND 
rC-PARP    WIT"    RFOUIPtC 
FLOW    RATE.        4N 
A'lXIL IA»Y     INLC*     MUST 
ftr     SUED     T-      4CCEPT 
ANY    FlCESS 
RCQUIREn   FLrw. 


I        OtN  =        I 
I  A!*vRMAX(?)*vn(?t  | 


t 
t 

AS  SUM  I  M*;  THE  AUX. 

INLET?  SLL?W  FLOW  Tf 

ENTER  SEFCRF  THE 

C1FPUSFP, 

CALCULATE     LOSSES     ANO 

TOTAL     PRESSURE    OF    THE 

r-mntNEn   FLOW. 


1       oc   = 

0.5"0(21          1 

1       1AIJ« 

'   or 

-    OIN        1 

cr   -  qaux      i 


ALCULATE     STATIC 
6ESSURE     I^MECIATELT 
FT.     OF    THE     L  IP. 


V12    •    OIK/AI 

Vfl     »    V!2/VC(2> 

PRL2    - 
TABLEtVRT.PRLT, 

ve.  jli 


I    O.S»RH->W*VI2*Vi,2       t 


alculate  cp^bineo 
low  pressures  ano 
ijx.   1ki.et  arfa. 
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)          auux 
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1          '/Mux 

=    0. 

1              °T<MjX 

-    0. 

BRftUX    =    0.^0 
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pb  »ux*on(  :>  ♦   spi 

HYP     *      PTVlX    -      5PI 

VIAU*     = 
*CRT(  ?.«0YO/PH]W( 


QAUX/V  MUX 


t 
I 

the  totu  phes^ure 
thf  co^in^o  flcw  i 
calculated  as  the 

••ASS    4E  IGmtpc    AVERA 
OF     THE    C CHAINING 
CLCMS. 


12            1       N"TE 
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1 
1 
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06 

CALCULATE    'HE 
INTERNAL    LENGTHS    TJ 
THE    DIFFU^''    ETPY. 


ELEHT    »    flcXT/o. 
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0!  I/ELTXTI 
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1.27374-MAUX" 
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04.12 >l 
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20 
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13  . 
112 


FLO  *  ELENT  . 
El  ♦  J.54401«02 
FLAUX 

ELI  •  S.5«01 

SIN  »  ELO 


«  ELL  .Gc.  ELO 


1 

ELN    =    ELL               1 

1 

1                       25 

I  / 

/  4.01 
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I 
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CALCULATE    ThE    IIP 
LOSSES    Fn«    FACH 
SITUATION. 
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I       OLIPIll    ■    1.    - 

I        TA8LEI WPT.PBLT, 
I  VELRIll.JLI 


1 

NUTE 

16 

. 

BEGIN 

00 

LOOP 

• 

l«    J    ■ 

3. 

1 

l< 

1 

NUN  8 

17 

0LIP1JI    *    1.    - 
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Appendix  E.   ASSUMED  CONSTANTS  AND  LIMITS 

It  was  necessary  to  make  numerous  assumptions  and 
impose  limits  on  various  quantities  in  the  subroutine.   Some 
of  the  assumptions  were  somewhat  arbitrary  and  may  be  mod- 
ified whenever  good  reasons  to  do  so  are  discovered. 

Some  of  those  assumptions  and  limits  are  discussed 
here  in  the  approximate  order  in  which  they  are  encountered 
in  SUBROUTINE  NACEL. 

The  first  assumption  that  should  be  mentioned 
is  that  everything  in  the  model  is  based  on  a  nacelle  with 
a  circular  cross  section.   No  attempt  was  made  to  account 
for  other  shapes. 

The  decimal  constant,  ZK,  which  determines  the 
amount  of  available  inlet  area  in  the  auxiliary  inlet  ring 
has  an  arbitrary  value  of  0.8.   The  real  case  may  vary 
widely  from  this  depending  on  the  details  of  the  design. 

The  incipient  cavitation  number  on  the  turning 
vanes,  SIGTV,  has  a  value  of  0.4.   This  value  was  accepted 
as  a  reasonable  one  for  the  speed  range  of  interest.   No 
examples  with  which  to  compare  were  found  in  the  literature. 

The  angle  of  attack  at  cruise,  TRIM(l),  has  been 

artificially  limited  to  a  maximum  of  three  degrees  because 

the  external  pressure  coefficients  are  not  single  valued 

with  respect  to  L/Dra  or  V±/Vm   for  c*  *  4°.   The  highest 

tabulated  value  for  angle  of  attack  is  6°  and  extrapolation 

would  probably  not  be  valid. 
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The  range  of  tabulated  velocity  ratios  is  from 
0.7  to  1.25.   Extrapolation  is  allowed  in  the  subroutine 
but  the  results  could  contain  large  errors  if  extrapolated 
very  far. 

The  largest  value  of  L/Dm  in  the  tables  is  2.50. 
This  is  also  allowed  to  be  extrapolated  and  the  same  remarks 
apply  as  for  the  velocity  ratios. 

D^/Dm  is  limited  to  0.90,  which  assumes  a  wall 
thickness  allowance  of  0.05Dm.   This  value  is  arbitrary  and 
may  be  modified  as  desired. 

The  pressure  loss  factor  for  the  auxiliary  inlet 
is  assumed  to  be  0.2  which  contains  an  underlying  assumption 
that  its  inlet  velocity  ratio  is  always  near  1.0. 

Diffuser  lengths  are  limited  to  those  correspond- 
ing to  double  angles  of  6°  to  20°.   See  chapter  4  for  an 
explanation  of  the  limits. 

In  the  wetted  surface  calculation,  the  area  of 
fairings  into  the  foils  and  strut,  as  well  as  the  missing 
area  in  way  of  the  foils  and  strut  are  ignored.  The  nacelle 
is  assumed  to  be  a  solid  body  of  revolution. 

Interference  drag  between  the  nacelle  and  the 
other  components  is  assumed  to  remain  constant  even  though 
the  nacelle  size  changes.  No  attempt  is  made  to  account 
for  the  effect  of  V./V  or  the  depth  of  submergence  on  drag. 
Chapter  8  contains  recommendations  concerning  this  subject. 
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The  weight  estimate,  as  discussed  in  chapter  6, 
is  only  a  rough  approximation  with  a  series  of  assumptions 
as  to  the  construction  and  material. 


-134- 


1347C 


"3 


Conner  ,..,„. 

A  study  of   ram  type 
ntake  parameters. 


5  SEP72 


DISPLAY 


Thesis 
C713 


Conner 

A  study  of    ram   tvpe 
intake   parameters. 


T34733 


thesC713 

A  study  of  ram  type  intake  parameters 


illliliii  llin  ll'.il  "I" 

3  2768  002  09328  8 

DUDLEY  KNOX  LIBRARY 


